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The interplay between intrinsic and surface/interface- induced magnetic anisotropies strongly in- 
fluences magnetization processes in nanomagnetic systems. We develop a micromagnetic theory to 
describe the field-driven reorientation in nanomagnets with cubic and uniaxial anisotropies. Spin 
configurations in competing phases and parameters of accompanying multidomain states are cal- 
culated as functions of the applied field and the magnetic anisotropies. The constructed magnetic 
phase diagrams allow to classify different types of the magnetization reversal and to provide de- 
tailed analysis of the switching processes in magnetic nanostructures. The calculated magnetization 
profiles of isolated domain walls show that the equilibrium parameters of such walls are extremely 
sensitive to applied magnetic field and values of the competing anisotropies and can vary in a broad 
range. For nanolayers with perpendicular anisotropy the geometrical parameters of stripe domains 
have been calculated as functions of a bias field. The results are applied to analyse the magnetiza- 
tion processes as observed in various nanosystems with competing anisotropies, mainly, in diluted 
magnetic semiconductor films (Ga,Mn)As. 
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I. INTRODUCTION 

In magnetic nanostructures complex physical processes 
on surfaces and interfaces give rise to enhanced uni- 
axial magnetic anisotropies [H [5J [31 3]. The inter- 
play between these induced and intrinsic (magnetocrys- 
talline) anisotropies strongly influences the magnetiza- 
tion processes in many important classes of nanoscaled 
systems, such as ferromagnet-antiferromagnet bilayers 
[3 H], thin epilayers of diluted magnetic semicon- 
ductors [9l Unj [TTJ [121 1131 IH! or in magnetic nanopar- 
ticles [ini Uni [13 [m U^i and is the reason of various 
remarkable effects involving complex spin reorientation 
[ini [3 [ini [HI 111 [13 [Ml [23 [23 [2Z1 [21] and the evolu- 
tion of specific multidomain states j29l [30l [311 [32l [33] . 

Most theoretical studies in this field are restricted to 
the investigations of Stoner- Wohlfarth processes through 
coherent switching [34 in models with uniaxial and cubic 
anisotropies 35, 36, 37, 38, 39, 40, 41J. Such theories de- 
scribe magnetization reversal in a limiting case of ideally 
hard magnetic materials. However, in real magnetic ma- 
terials the reversal processes will usually take place by 
the formation of heterogeneous states consisting of the 
competing phases and their transformation under the in- 
fiuence of the applied field. The analysis of such mul- 
tidomain states in systems with competing anisotropies 
and their influence on magnetization reversal is the sub- 
ject of this paper. These investigations can be executed 
within a regular micromagnetic theory [42' '43] , adapted 
to nanoscaled systems (see e. g. Refs.l44, 45, 46). 

In section II we introduce the phenomenological model 
and methods; in the next section we derive all pos- 
sible magnetic configurations in the system, calculate 
their stability limits, and describe reorientation effects 
[m [37] HE]- These results enable us to calculate 



the parameters of the multidomain states and analyse 
the magnetization processes (section IV). In section V 
we apply our results to interpret reorientation effects 
and magnetization reversal as observed in experimental 
works on nanolayers of diluted magnetic semiconductors, 
FM/AFM bilayers, thin films of Heusler alloys, and mag- 
netic nanoparticles. In section VI the calculated equilib- 
rium parameters of the isolated domain walls and stripe 
domains are used to analyse recent experimental results 
in (Ga,Mn)As films with perpendicular anisotropy. 

II. PHENOMENOLOGICAL MODEL 

Within the standard phenomenological theory the 
magnetic energy of a nanoscale ferromagnetic sample can 
be written as a functional Wm — J w(r)dV with an en- 
ergy density 

^-^E(|^)'-M-H(^) - iM.HW-fz.,, (1) 

where m = M/Mq (Mg = |M|) is the normalized value 
of the magnetization vector M, A is the exchange con- 
stant, and H^'^^ and H'^'^^ are the external and demagne- 
tizing fields, respectively. The anisotropy energy density 
includes uniaxial anisotropy (K^) with the axis a and 
cubic anisotropy {Kc) with unity vectors nj along cubic 
axes 

1 ^ 

u-a(M) = -X„(m.a)'--if,^(m.nj)4. (2) 

The coefficients and Kc are assumed to be positive. 
Hence, a and iij directions are easy uniaxial and cubic 
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magnetizaton axes, respectively. 

The equilibrium distribution of the magnetization 
M(r) is generally spatially inhomogeneous. It can be de- 
rived directly by solving the equations minimizing the en- 
ergy functional Eq. (jlj together with the Maxwell equa- 
tions. Thus, the micromagnetic problem is formulated as 
a set of non-linear integro-differential equations [43 . In 
many classes of magnetic systems a strongly pronounced 
hierarchy of magnetic interaction scales allows to reduce 
the micromagnetic problem to a set of auxiliary simpli- 
fied problems |15J |33]. The procedure includes: (i) the 
calculation of spatially homogeneous equilibrium states 
by minimizing energy 

u;o(M) = -M • H + Wa{M) (3) 

in an (internal) magnetic field H for fixed values of the 
material parameters in Eq. (|3|. The solutions of (i) are 
used to construct magnetic phase diagrams in compo- 
nents of the external magnetic field (ii) and to calculate 
the equilibrium parameters of multidomain patterns and 
the structure of domain walls (iii). 

In the rest of the paper, we apply this program to 
the model given by Eqs. ([T]) and In order to make 
transparent the representation of internal homogeneous 
states and the phase diagrams, we restrict our discussion 
to the case of co-planar arrangements of easy axes and 
applied fields. Generalizations of this model are discussed 
to the end of the next section. 



III. REORIENTATION TRANSITIONS AND 
METASTABLE STATES 

In many cases of practical interest, the direction of the 
uniaxial anisotropy a lies in the plane spanned by two of 
the cubic axes nj (see Eq. To be specific we define 

this plane as xOz plane assuming that z is directed along 
a. In this case energy wq from Eq. ^ can be written as 
a function of the angle 9 between M and a. Introducing 
the reduced energy ^{9) = ■wq/Kc + 1/8 one obtains 

^(9) = -- cos 4(61 -a)- >c cos^ 9 - h cos{9 - iP) , (4) 
8 

where 

x=K^IK,, h = H/H,, Hc = K,Ma, (5) 

a is the angle between the unaxial a and cubic rii axis, 
the angle tp defines the deviation of the magnetic field H 
from the easy axis a in the xOz-plane; correspondingly, 
hz = h cos ip is the reduced field component along the 
uniaxial easy direction, and — hsinip is the perpen- 
dicular component. 

Energy Q is a function of variable 9 and includes four 
material (control) parameters, namely, angle a, ratio >c 
and reduced magnetic field components, h^, h^. Model 
Q has been introduced in 1964 by Torok et al. [3S] 
for ferromagnetic films with misorientated uniaxial and 



biaxial easy magnetization directions. Previous inves- 
tigations of Q have been restricted to limiting cases 
of a = and a ^ n/A and were mostly concentrated 
on investigations of coherent rotation processes (Stoner- 
Wohlfarth regime) (see Ref. i34jand bibliography in Ref. 
l43|l . Within this approach switching processes are identi- 
fied with the boundaries of the metastable states (critical 
astroids). In this section we give a comprehensive analy- 
sis of model ^ in the full range of the control parameters 
{a, x^hjtp). In particular, we show that the analysis of 
the metastable states only is not sufficient for the under- 
standing of magnetization reversal in nanosystems with 
competing anisotropies. The peculiar evolution of the po- 
tential profile ^ under the influence of the applied field 
and specific reorientation effects are found to be crucial 
for the magnetization processes in this class of magnetic 
materials. 

The stationary solutions with the equilibrium values of 
9 are derived from the equation <i>e = 0: 

/isin(6' - V) = -^sin4(6' - a) - >fsin26' (6) 

(Here we introduce a common notation for derivatives 
f{x)xxk = d^f/dx^). The equation for the lability lines 
of the solutions, = 0, reads 

hcos{9 - V") = -2 cos 4(6* - a) - 2>^cos26' (7) 

and determines the stability boundaries of the solutions 
together with Eq. (|6| . Critical points of the transitions 
are determined by the set of equations <i>exfc = , fc = 
1, 2, 3. The degeneracy of the solutions, 9i with i — 1...L, 
that provide global energy minima in the system, 

H9i)^H92) = ... = H9l), (8) 

determines the coexistence regions between L magnetic 
phases and the conditions for first-order transitions in the 
magnetic phase diagrams. It is convenient to present the 
solutions 9{hx, hz, X, a) and critical regions for Q in a 
set of phase planes {hx,hz) parametrized by the factors 
a and x. We start our analysis from the limiting cases 
a = and 7r/4, and after that highlight main features of 
the general model 

For symmetric cases with the easy-axes orientations 
along one the cubic axis (a = 0) and along the diagonals 
between them (a = 7r/4) the Eq. Q can be written as 

$(6») = ipi cos(46i) - >c cos^ 9~h cos{9 - V') , (9) 
8 

with "— " for a = and for a = 7r/4. Gener- 
ally model ([9]) describes magnetic states in a planar fer- 
romagnet with competing uniaxial (second-order) and 
biaxial (fourth-order) magnetic anisotropy. The model 
has been applied for many bulk and nanoscale magnetic 
systems, including reorientation effects in rare-earth or- 
thoferrites j42j |49] , several classes of intermetallic com- 
pounds |50j . first-order magnetization processes in high- 
anisotropy materials |38j . and for magnetic nanolay- 
ers with surface/interface-induced magnetic anisotropy 
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[1 H la IMl SOI SSI iZ] ■ The model from Eq. ^ has also 
proved to be valid for diluted magnetic semiconductors 
as a novel class of magnetic materials El UHl ITTl I48j . 

The invariance of the potential Eq. (l9| under the trans- 
formation 

x^-x, + + 71-/2, (10) 

means that the cases with different sign of x trans- 
form into each other by rotation of the reference system 
through an angle 7r/2. This invariance allows one to re- 
strict the analysis to positive values of x. The analysis 
of (|9| yields four topologically different types of (/i^, h^)- 
phase diagrams depending on values >f > 0, as shown in 
Fig-E 



Under transformation (10 1 the equations of equilib- 
rium (|6|, (0 for the potential ^ with a = are con- 
verted into those equations for a = 7r/4. Thus, for the 
same values of x the lability lines for both cases trans- 
form into each other by a rotation through 7r/2 (Fig. 1). 

For X > 5 the lability lines have a similar shape as 
the Stoner- Wohlfarth astroid [34 . In the limit of large 
X, the lability line asymptotically coincides with this as- 
troid for simple uniaxial ferromagnets. As x decreases 
from 5 to zero the lability lines transform into an eight- 
cusp line with the shape of the classical wind rose. At 
the parameter value x = 5, a, bifurcation takes place 
by the formation of so-called "swallow tails" in one 
pair of opposite cusps (Fig. 1(B), (F)). In the interval 
5 > X > 1, the swallow tails gradually widen, and at 
X = 1 the lower cusp points reach the horizontal axis. In 
the interval 1 > x > the phase diagrams have regions 
with overlapping swallow tails (Fig. 1(C), (G)). Finally 
at X = 0, that is for zero uniaxial anisotropy, the lability 
line transforms into the wind rose with eight corners. For 
the two special orientations of the uniaxial anisotropy 
in Eq. (|9| with a = and a = n/i the lability lines 
are identical. Still, the phase diagrams are fundamen- 
tally different as they pertain to different magnetic states 
and different coexistence regions of metastable magnetic 
states. Depending on the control parameters, there are 
regions with L = 2,3, or 4 degenerate states, and, con- 
sequently, first-order transitions involving two, three and 
four phases. 

For X > 5 the first-order lines between two magnetic 
states are segments of straight lines connecting opposite 
cusp points h, - 0, {\h,\ < \h'^\ = 2(xT 1) (Fig. [ijA), 
(E)) [46 . At zero field the transitions occur between 
antiparallel magnetic states 9i — and 02 — tt. For 
finite transversal magnetic fields \hx\ < h% the solutions 
for coexisting phases are determined from Eq. (|6| with 
hz^Q and a = (7r/4), 



s n^ e s n 61 ± — 

2 4 



0. 



(11) 



These solutions describe a gradual decrease of the mag- 
netization component m^. In the endpoints of the first- 
order transitions = ±ft,J^ the magnetization vectors in 
both phases are perpendicular to the easy direction. 



For 5 > X > the evolutions of the magnetic states 
within the swallow tails are different for the two models 
(see potential profiles in Fig. [2]). For the model with 
a — 0, the potential wells corresponding to the global 
energy minima are swapped within the swallow tails (Fig. 
|2jA)). Hence, different canted states become degenerate 
in equilibrium along lines of first-order transitions (lines 
fli^i, ai62, 0263, 0264 in Fig. |2jA)). These lines meet 
the transition line aia2 between symmetric phases, 9i, 
62 = — 01, ya the points ai and 02 where three phases 
coexist. The coordinates of points ai and 02 are [50] 



h,j; ^ ±2 sin {cos 29 + x), ^ 
and the solutions for the coexisting phases read 



h=±0, 02 = TTTi 



±7r/2, 



(12) 



(13) 



where sm9= {-1 + Vl + 3>f)/3. 

For the model with a = 7r/4, the minimum in the po- 
tential is unique for the swallow tails in the parameter 
range 5 > x > 1. Hence, the transformation of the en- 
ergy profile involves only metastable states (Fig. |2][B)). 
However, the particular transformation of the metastable 
states in these regions plays an important role in the 
evolution of the domain wall profiles (see Sec. VI). 
In the interval I > x > the first order transitions 
arise within the region of the overlapping swallow tails. 
The transition line is a segment C1C2, with the points 
ci = (0;2(x-l)), C2 = (0;-2(x-l)). Along this 
line segment two phases coexist with solutions 0i and 
02 = —6*1. The solutions 9i are given by the equation 
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cos 



0, 



(14) 



that can be derived from Eq. ([6|. The first-order tran- 
sition line from ci to C2 crosses the other transition line 
along (/ij.; 0) in the origin. Hence, in this point four mag- 
netic phases coexist with 
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■ arccos x, 02 



(15) 



For X = (zero unaxial anisotropy) both potentials 
([9]) are converted into the model of a cubic ferromagnet. 
The corresponding phase diagrams (Fig. [TJD) and (H)) 
are identical and include two lines of first-order phase 
transitions between symmetric states. In the origin four 
degenerate phases with magnetization along the cubic 
easy axes rii and 112 coexist. 

In the general case with a misalignment between uni- 
axial and cubic easy axes given by the parameter a the 
potential $ (Eq. Q) is a periodic function of a with peri- 
odicity 7r/2. Thus an analysis in the range < a < tt/A 
covers all physically different states. Here we describe 
the evolution of the (h^, hz) diagrams when a varies from 
zero to 7r/4. The sets of diagrams in Figs. |3] and |4] show 
the transformation of the transition and lability lines. 
The case with nonoverlapping swallow tails for the pa- 
rameter range 5 > x > 1 is presented in Fig. |3] For 
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FIG. 1: (Color online) The phase diagrams of magnetic states in components of the internal magnetic field {hx, /iz) and 
different values of the parameter >c for the ratio of uniaxial and cubic anisotropy (Eq[5|. The upper panel (A)-(D) is for 
systems with a = Q, the bottom panel (E)-(H) is for systems with a — 7r/4. Two-headed vectors show orientations of the 
uniaxial axis a and the cubic axes. The plots present the topologically different types of phase diagrams: (A), (E) x > 5, 
(B), (F) 5 > x: > 1, (C), (G) 1 > >c > 0, (D), (H) x = 0. Black lines indicate stability limits of metastable states. Red lines 
give the first-order transitions between different magnetic phases (see text for details). 
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FIG. 2: The {h^, h^) phase diagrams for }< — 0.7: a — (A), a — n/A (B). Potential profiles ^{0) are sketched for various 
points in the phase diagrams to illustrate the evolution of the magnetic states in both models. 



small a{>c) two lines of the phase transitions between 
canted phases still exist (Fig. [sj^B)). With increasing a 
the points ai and 02 for the three-phase coexistence move 
to either of the cusp points 62 and &3 (Fig. [3][C)) Af- 
ter these points have merged only two-phase transition 
lines exist in the system (line 6164 in Fig. [3]jD), (E)). 
Thin (blue) lines in Fig. |3|B)-(F) indicate the values of 
the magnetic fields where two metastable states have the 
same energy. They are not connected with any physi- 
cal processes in the system but help to understand the 
transformation of the energy profiles. 

The phase diagram with overlapping swallow tails for 
the parameter range 1 > >f > is shown in Fig.[4j In this 
case the transition lines between pairs of canted phases 



in the limit a — 0, Fig. |4][A), gradually transform into 
straight line segments for transitions between pairs of the 
symmetric phases in the limit a = 7r/4, Fig.|4][D). During 
this process the points of the three-phase coexistence ai 
and a2 move towards each other (Fig. QB), (C)), and 
merge into the point with four-phase coexistence at the 
origin for a = n/4, Fig. |4]^D). 

The sets of modified astroids in Figs. [T] and |3] repre- 
sent geometrical singularities studied by a special field 
in mathematics known as catastrophe theory |51j . It was 
found that for rather general form of potentials there ex- 
ist only seven fundamental types of singularities referred 
to as catastrophes |51j . Four of them are realized in the 
lability lines of Figs. [l] and [3] The astroid lines, where one 
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FIG. 3: {hx, hz) phase diagrams for x — 1.2 and different 
values of a demonstrate the transformation between the two 
symmetric cases with a = (A) and a = 7r/4 (F). 



in micromagnetism. 

The transformation of the common Stoner-Wohlfarth 
astroid (as 4-cusped hypercycloid- Fig. [ijA), (E)), into 
the 8-cusped hypercycloid with the wind rose shape (Fig. 
[IJD), (H)) occurs in many magnetic systems with com- 
peting anisotropics and has been investigated during the 
last forty years. To the best of our knowledge the 8- 
cusped hypercycloid has been firstly obtained in Refs. 
123 [S2] (see also the remarks about earlier conference 
contributions in Ref. I35B. The transformation from the 
common astroids (Fig. IT]) into a lability line with swal- 
low tails, and the further evolution of these curves to 
the wind rose has been obtained in Ref. EHl Torok et 
al. also demonstrated several diagrams of lability line 
for model with misorientated anisotropy axes (Eq.([9])). 
In many following papers (see, e.g. Refs. [3S1 133 TmI 
pecularities of lability lines for the model (|9| have been 
investigated. The coordinates of the critical points for 
([9]) were calculated in Refs. 50, 54, 55, 56 The solutions 
for the first-order phase transition lines and the coexist- 
ing states have been carried out in Refs. [571 HHI (see also 
Ref. Sni)- In this paper we have given an exhaustive 
summary of model 

IV. MAGNETIC PHASE DIAGRAMS 
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FIG. 4: {hx, hz) phase diagrams similar to that in the pre- 
vious figure but for >c = 0.7. 



local minimum merges with a local maximum, are named 
fold catastrophes. The edge points where two folds meet 
have two minima merging with a maximum. These sin- 
gularities are known as cusp catastrophes. By a charac- 
teristic discord in the terminology, the feature known in 
magnetics as "swallow tails" as shown in Fig. 1(B), (F) 
are called butterfly catastrophes in mathematics, while 
the triangular regions, as those with the cusps €1,62 in 
Fig. [S^D), are called swallow tail catastrophes. The la- 
bility lines in Fig. [T] belong to a family of hypercycloids. 
In present article we shall adhere to the terminology used 



In the previous section the solutions for possible mag- 
netic states have been presented as functions of the in- 
ternal field. For ellipsoidal magnets with a homogeneous 
magnetization M(h) the equation [43] 

h^*^' = h + 47rX^"iNm(h) (16) 

establishes the correspondence between magnetic phase 
diagrams in terms of the internal field h and those in 
terms of the external field h'"^) (N is the demagnetizing 
tensor). For phase diagrams in Figs. [l]and|3]the phase 
diagrams in external magnetic field components are plot- 
ted in Figs. [5] and [6] correspondingly. 

Thin lines in Figs. [5] [6] define values of the external 
fields in which the internal field within the stable phases 
reaches the boundaries of the metastable region. The 
transition lines (hf^) in Figs, [l] |3] transform into the ar- 
eas bounded by thick (red) lines in Figs. |5][6j These areas 
define maximum possible regions where thermodynami- 
cally stable multidomain states of the competing phases 
can exist ji^ 1151 155]. 

For x: > 5 and a = the {h^, h^) diagram in Fig. [l|A) 
converts into that in Fig. |5jA). The transition line in 
Fig.jljA) transforms into an area of a two-phase multido- 
main state that is bounded by an ellipse (red line). For 
I >«:| < 5 and a = the diagrams in the terms of external- 
field components become rather complicated: Fig.JSiB) 
is obtained by mapping the phase diagram in Fig. [iTB), 
and Fig. [5](C) by mapping Fig. [iJD). The lines of the 
two-phase transitions in Fig. [iJB), (C) convert into ar- 
eas of two-phase domain structures (DS); the points of 
three-phase coexistence (fT2| "swell" into the triangular 
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FIG. 5: Magnetic phase diagrams in the components of the external field hx , hi for a = ( (A)- (C)) and a = 7r/4 ((D)- 
(F)). Thick lines limit regions of three- and four-phase (blue) and two-phase (red) multidomain states. Arrows show magnetic 
configurations in the (co-existing) domains. All calculations have been carried out for a spherical sample {Na = 1/3). 
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FIG. 6: Magnetic phase diagrams in the components of the external field hx , hi for x — 1.2: (A)- a = 0; (B)- a = n/15; 
(C)- a = tt/8; (D)- a = 7r/4 and x = 0.7: (E)- a = 0; (F)- a = 7r/15; (G)- a = 7r/8; (H)- a = tv/4. Thin dotted lines are the 
lines of the constant internal field and the constant phase fractions. 



regions of three-phase domain structure; and the point 
(0,0) in Fig. [iJD), where four phase 6j — i^j/^ coex- 
ist, transforms into a rectangular area with a four-phase 
multidomain state. 

The phase diagrams in Fig.[l]jF), (G), (H) for systems 
with a = 7r/4 are mapped into the phase diagrams in 



Fig. [sj^D), (E), (F), correspondingly. For x > 1 the 
phase diagram in Fig. [sJD) includes one region of two- 
phase multidomain states and two swallowtail pockets 
in the metastable region. For 1 > xr > h^'^^-phase 
diagrams include a rectangular area, where four-phase 
domain states are stable with spin configurations in the 
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domains described by Eq. (15 1. Adjacent to this area, 



there are four regions in the phase diagrams with two- 
phase multidomain states (Fig. |5[E)). Finally the h^*^^- 
phase diagram for >f in Fig.^F) becomes identical 
to that in Fig. Isj^ 



For a varying from zero to 7r/4 the phase diagrams are 
plotted in the case k=\.2 (Fig. |6];A)-(D)) and h = 0.7 
(Fig. [6]jE)-(H)) and reflect the complex transformation 
of the regions of the DS existence. For yc ~ 1.2 the 
triangular areas of three phase DS and the regions of two 
phase DS deform (Fig. |6jB)) and then disappear at all 
(Fig. [6][C)) leading to a phase diagram with one distorted 
ellipse of two phase DS. The phase diagram for h — 0.7 
demonstrate another kind of transformation. Now, the 
regions of three-phase DS do not disappear (Fig. [6](F), 
(G)) but, on the contrary, join to form a rectangular area 
of four phase DS (Fig. [6][H)). Here in Fig. [6] the thin 
dashed lines denote the lines of constant internal field and 
constant phase fractions (in regions of stable two-phase 
DS). When an external magnetic field is varied following 
these lines, then either domain walls are displaced or the 
magnetization rotates in each domain, respectively. 

In the next section we apply the diagrams of solu- 
tions in Figs. 1|2|3|4 and the phase diagrams in Figs. 
|5|6| to analyze the magnetization processes in nanosys- 
tems with competing anisotropics. We emphasize that 
those diagrams describe two limiting cases of ideally hard 
and ideally soft magnetic behavior. In ideally hard mag- 
netic materials magnetization processes occur via evolu- 
tion of metastable states. Magnetic phases exist every- 
where in their stability regions up to their boundaries 
{Stoner-Wohlfarth regime). The corresponding magneti- 
zation curves (dotted lines in Fig. [t]) are characterized 
by the widest possible hysteresis cycles [13], and single 
domain states are realized in these systems. In the op- 
posite case of ideally soft magnetic materials the mag- 
netization reversal occurs via the evolution of thermo- 
dynamically stable states. Such anhysteretic processes 
involve the formation of multidomain patterns. These 
spatially inhomogcncous states are composed of domains 
formed by the competing phases of the magnetic-field in- 
duced first-order transition [42. 43^ . Extended regions of 
multidomain states have been observed during reorien- 
tation processes in several groups of bulk magnetic sys- 
tems (e. g. orthoferrites and easy-axis antiferromagnets 
P^IBU] ). For these magnetically soft, low anisotropy sys- 
tems the multidomain states are described by the phase 
theory equations [42l |43] . The phase theory approxima- 
tion stricly is valid only if the characteristic sizes of the 
sample are much larger than the sizes of domains, and 
transitional regions between domains are localized into 
narrow domain walls [43| . Both these requiments are met 
only in ideally soft, massive magnetic samples. Thus, in 
soft magnetic materials the magnetization processes are 
mainly determined by occurrence of the first-order phase 
transitions and the evolution of the magnetic states in the 
coexistence phases during these transitions (solid lines 
m Fig. [t]). The magnetization processes in real mate- 



rials are intermediate between these two limiting cases 
and include both evolution of the metastable and mul- 
tidomain states. In magnetic nanolayers domain sizes 
usually exceed the layer thickness. In magnetic nanopar- 
ticles only few domain walls are observed, and in suffi- 
ciently small particles multidomain states are completely 
blocked. On the other hand, coercivity of the magnetic 
nanosystems prevents the formation of the equilibrium 
states and causes hysteretic magnetization reversal. 



V. COMPARISON WITH EXPERIMENT 



The phase diagrams of solutions in Figs. [T]|4]can be ap- 
plied for explanation of magnetization processes in many 
nanomagnetic systems with competing anisotropics, for 
example, in thin films of diluted magnetic semiconductors 
(DMS), in ferromagnetic(FM)/antiferromagnetic(AFM) 
bilayers [5J[23], Heusler alloys pT|, and/or nanoparticles 
[13117]. First, we consider the phase diagrams with sym- 
metric arrangement of easy uniaxial and cubic anisotropy 
axes (Fig. lj2l, as applied for nanolayers of DMS. And 
then we give examples of systems with different values of 
angle a between anisotropy axes. 

Layers of diluted magnetic semiconductors represent a 
new class of materials with a strongly pronounced com- 
peting character of the magnetic anisotropy. In exist- 
ing (Ga, Mn) As nanolayers the ratio x of uniaxial and 
cubic anisotropy varies in a broad range and is con- 
trolled by strains, temperature and hole concentration 
[g Unj ifH Ha [Til [Ml |23]. The magnetization pro- 
cesses in Gai_a;Mna,As thin films grown by molecular 
beam epitaxy on GaAs(OOl) substrates are described by 
the diagrams of solutions for highly symmetric geometry, 
a = 0;7r/4 (Fig. [T|). The in-plane magnetization rever- 
sals in these systems are determined by the competition 
of cubic anisotropy with easy axes < 100 > and uniax- 
ial anisotropy favouring the directions of < 110 > type. 
Thus, the solutions of Q for a = 7r/4 are applicable in 
this case (Fig. [iJE)-(H)). The main features of in-plane 
magnetization processes in such layers are summarized 
in Fig. |2JB) and were experimentally investigated in a 
number of works [H [21 [13 [53] . 

In Ref. [29] the 300nm thick Gai_xMnxAs(x = 0.03) 
samples were studied combining direct imaging of mag- 
netic domains and SQUID magnetomctry. At tem- 
peratures above 30K the samples exhibit the uniaxial 
anisotropy with easy axis along [110], whereas for tem- 
peratures below 30K the magnetization vector deviates 
from this direction indicating the dominance of the four- 
fold symmetry. According to our phenomenology these 
magnetic films followed the temperature transition from 
the phase diagram in Fig. [iJF) for dominating uniaxial 
anisotropy (>f > 1) to that in Fig. [TJG) with competing 
character of anisotropy (>f < 1). The angle between the 
magnetization and the axis [110] is determined by the 
Eq. (fl5|. For T=15K (the ratio x = 0.42 was obtained 



from fits of the hard-axis magnetization curves and cal- 
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culating the easy-axis orientation from Eq.([9|)) this an- 
gle is 32° which agrees with experimental results. The 
magnetization processes for high temperature (>f > 1, 
Fig. [iJF), [sj^D)) proceed through the nucleation and ex- 
pansion of domains with two orientations of the mag- 
netization vector. In fields applied along the easy axis 
[110] the evolution of the DS is accompanied only by 
the 180° domain wall movement (Fig. [6|^D)), and the 
metastable states in swallow tails (see energy profiles in 
Fig. [2]jB)) can be considered as the nuclei of domains. 
During the magnetization processes along [100] axis the 
domains not only nucleate and expand but the magneti- 
zation rotates inside each domain of various phases (Fig. 
p[D)). The magnetization reversal for low temperature 
(>f = 0.42) along one of the cubic easy axes proceeds 
in three stages [29 through the formation of intermedi- 
ate domains (Fig. [l|^G), [sj^E)). In the first (and the 
last) stage a transversely magnetized domain nucleates 
indicating the entering into the area of two-phase DS 
in Fig. [sjjE). With field increasing the completely re- 
versed domains nucleate and propagate rapidly through 
the sample indicating the beginning of the area with four 
phase DS in Fig. |5]^E) ^5]. It is remarkable that the one 
stage switching processes are also possible and are ac- 
companied by the transformation of four phase domain 
structure for some directions of magnetic field (blue open 
points in Fig. [sJE)). The magnetization curves in Fig. 
[7](A), (B) are typical for the in-plane geometry display- 
ing biaxial character of the anisotropy. The successive 
switching of the magnetization in Fig. [7]^A) are caused 
by the redistribution of the metastable minima in the en- 
ergy profiles for varying magnetic field (Fig. [2]^B)). The 
hysteresis loops of such a type are more pronounced for 
purely cubic anisotropy {>c = 0) and were observed for 
603 nm-Gao.957Mno.o43As films [50]. If the applied mag- 
netic field makes the angle with the [100] axis in the range 
(0; 7r/4), and initially spins are aligned along [100], then 
the first incoherent reversal is related to the appearance 
of domains with [010] magnetization, whereas the second 
jump IS due to the [100] domain (Fig. [TJD)). 

In Ref. [SU the character of in-plane magnetic 
anisotropy has been determined by means of transport 
measurements. All layers were patterned into 40-60/im 
wide Hall bar structures, and a strong anisotropic mag- 
netoresistance effect provides a very convenient method 
to study the anisotropy at fixed temperature. The resis- 
tance polar plots of transport measurements for preva- 
lent biaxial anisotropy [55] look similar to the phase di- 
agram in Fig. [sJC). The [110] uniaxial anisotropy leads 
to the narrowing and subsequent disappearance of the 
four phase DS area (Fig. [5]^D), (E), (F)). As well, it was 
shown that an additional uniaxial anisotropy with [010] 
easy axis is present in the system. This anisotropy re- 
sults in the formation of a two phase DS region splitting 
the rectangle with four phases (Figs. [sj^B), (C)). 

The solutions with a = (Fig. [i];A)-(D)). are realized 
for out-of-plane magnetic field and in-plane orientation 
of the magnetization [5TJ [571 HH]- The lability lines of 



phase diagrams for a = (Figs. [ijB), (F)) are similar 
to those for a = 7r/4 but the magnetization processes are 
quite different. The triple point with three coexisting 
phases inside the swallow tail (Fig. [TJB)) has a crucial 
influence on the magnetization reversal and is the reason 
of specific double shifted hysteresis loops (Fig. |7][C)) 
observed in many works [T^l |57j . 

In Ref. 127i Gai_xMnxAs films grown on hybrid 
ZnSe/GaAs substrates with a low Mn concentration {x w 
0.01) were chosen to identify the role of both types of 
anisotropies in the magnetic reversal process. Varying 
the hole concentration p and temperature T the ratios 
XT according to Figs. [l|A)-(D) can be swept. For the 
hole concentration p = 8.5 • lO^^cm""^ the temperature 
progression results in the succession of phase diagrams, 
namely. Fig. [i|A) for T=20K, Fig. [TJB) for T=7K, Fig. 
[TJC) for T=3K and Fig. QD) for T=1.5K. The switch- 
ing processes for high temperature exhibit the typical be- 
havior of a specimen magnetized along the hard direction 
(Fig. [TJA)). In this case, a domain structure exists with 
magnetization vector tilted with respect to the magnetic 
field. As the temperature is lowered the triple point in 
Fig. [iJB) denotes the existence of an additional stable 
magnetization state along the magnetic field direction. 
Thus, the subloops of the hysteresis curves (Fig. |7][C)) 
refiect the jump of the magnetization into this minimum 
accompanied by the three phase domain structure. The 
variation of magnetic field in the region spanned by the 
swallow tails lead to various scenarios of the DS transfor- 
mation (see energy profiles in Fig. |2]^A)). In particular, 
different cases (shown in Fig. [sj^B) by the red and blue 
open points) of the transition from multidomain states 
into a single domain state can be realized. For some di- 
rections of the magnetic field crossing regions with two- 
and three-phase DS one can observe even more com- 
plex hysteresis loops consisting of three subloops (Fig. 
|7|^D)). In our phenomenology the maximum hysteresis 
loops encircle the anhysteretic magnetization curves with 
three and two phase DS. Experimentally that kind of 
magnetization processes was observed in Co2MnSi and 
Co2MnGe Heusler alloys pT. With the temperature de- 
creasing the subloops in Fig. [tJ^C) broaden around the 
two steps of the magnetization, and a hysteresis loop with 
only a weak double shift is observed (at T=1.5K, the ex- 
perimentally measured value of k « 0.26 corresponds 
to phase diagram FigjTfE)). Note, that one should dis- 
tinguish the hysteresis loops of such type for the cases 
with a = 7r/4 (Fig. ^k)) and a = (Fig. [tJE)) 
because the magnetization processes are fundamentally 
different. For (Ga,Mn)As systems with low hole concen- 
tration {p — 3.0 • lO^^cm"^) uniaxial anisotropy is much 
larger than the cubic anisotropy |27]. This situation is 
described by model Q with K>b (Fig. [ijA), [sjA)). If 
the magnetic field is applied along the easy axis of uni- 
axial anisotropy (Fig. [ljB)j5]^B)) then the two phase DS 
transforms into a single phase state. As a remnant of 
the DS, 360° domain walls may remain in this state and 
can act as nuclei of new reverse domains when lowering 




or changing the external magnetic field. Experimentally 
such situations have been studied in magnetic field per- 
pendicular to the film and for out-of-plane magnetization 
vector [55] , 

In FM/AFM bilayers of cubic materials the intrin- 
sic cubic < 100 > anisotropy may compete with the 
uniaxial anisotropy induced by the exchange couplings 
between two layers To establish the exchange 

bias uniaxial anisotropy, the bilayer film is cooled in 
an in-plane magnetic field which determines the easy 
axis of induced anisotropy. In Ref. the magneti- 
zation reversal has been studied in an exchange-biased 
CaMnOs/Lao.erSro.asMnOs bilayer film grown on vici- 
nal SrTiOs < 100 > with the angle between cubic and 
uniaxial easy axes being a — 30°. With temperature de- 
creasing the magnetic films followed the transition from 
the phase diagram of solutions in Figs. [3|E) (T=160K) to 
that in Fig. QC) (T=5K). The magnetization processes 
for high temperatures involve only the redistribution of 
a two phase DS (Fig. [6|C)). The metastable states in- 
side each swallow tail lead only to the slight deformation 
of energy profiles and do not infiuence significantly the 
magnetization reversal (Fig. [3][E)). For low temperatures 
those metastable states become stable (Fig. |4]^C)) and 
alternatively a domain state can be realized. The three 
and two phase DS (Fig. |6jF),(G)) in the system result 
in complex hysteresis loops with a hint of a double shift 
(Fig. ^E)) [23. 

Double shifted magnetization curves with strongly pro- 
nounced subloops (Fig. |7|C)) and the astroid with swal- 
low tails (Fig. d^A)) were observed in Ref. 8 for metallic 
multilayersamples with the structure Si(100)/Cu(15nm)/ 
Ni8oFe2o(35nm)/NiMn(50nm)/Co(tnm)/Pd(15nm) 
grown by an e-beam evaporation system. The thickness 
was varied between 5 nm and 25 nm. It was shown that 
the double-shifted hysteresis loops (and parameter xr) 
could be tuned by several parameters, e.g., the variation 
of Co film thickness, and the field-annealed time. 

In magnetic single-domain nanoparticles the com- 
petition of the uniaxial anisotropy due to the en- 
hanced surface interactions and intrinsic magnetocrys- 
talline anisotropy stabilizes different multiple magnetic 
states in the system with the possibility of switching be- 
tween them (for details see Ref. [T5)) . Different geome- 
tries of relative easy axes alignment are realized in these 



nanoobjects. Phase diagrams with swallow tails (Fig. 
|6jC)) for misaligned easy anisotropy axes have been ob- 
tained for FeCuB nanoparticles [16] and Co clusters [17]. 
Astroids with rounded corners (Fig. [SjA)) have been ob- 
served in fcc-Co [ISj and BaFeCoTiO nanoparticles [19] . 

In general, in many cases of practical interest the geom- 
etry with noncoplanar competing anisotropics is observed 
[231 mi [53] . For example, in 50 nm thick Gao.91Mno.09As 
thin films grown on (311)A and (311)B substrates |23| the 
uniaxial anisotropics with [Oil] (or [233]) and [311] easy 
axes compete with cubic anisotropy of < 100 > type. In 
this case three dimensional phase diagrams, parametrized 
by the varios ratios of anisotropy coefficients and mag- 
netic field components, have to be constructed instead 
of 2D diagrams of solutions. Even for considered DMS 
films with (001) orientation, the out-of-plane magnetiza- 
tion processes can be considered as coplanar only with 
some restrictions. Indeed, for biaxial in-plane anisotropy 
and magnetic field perpendicular to the film one generally 
has a non-coplanar arrangement of the magnetization in 
domains. In that case, the phase diagram in Figj5](B) is 
only the cross-section of that more complex 3D phase di- 
agram. But due to the degeneracy of in-plane cubic [100] 
and [010] axes with respect to the magnetic field the main 
peculiarities of the switching processes can be readily ex- 
plained with the simple 2D phase diagram (Figl|B)). 
Therefore, the magnetic anisotropy geometry and mag- 
netic field orientation determine which phase diagrams 
of solutions (2D or 3D) is applicable in a particular case. 



VI. MULTIDOMAIN PATTERNS 

Multidomain patterns have been observed in a number 
of systems with in-plane |29j and out-of-plane magneti- 
zation ISHIMllSSllBllMllMlinS]- Particularly, iso- 
lated domain walls trapped in micropatterned constric- 
tions of (Ga,Mn)As films demonstrate a number of re- 
markable properties [67l [Ml IMl EQ] and can be used in 
differerent nanoelectronic devices (e. g. as magnetore- 
sistive elements) [57]. The fine structure of the isolated 
domain wall is of prime importance when different types 
of domain walls are observed |71j . Here in particular, we 
demonstrate that for the considered systems with com- 
peting anisotropics various types of domain walls exist 
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with large sensitivity of their appearance on material pa- 
rameters and external fields. Using the results of the 
two previous sections we calculate the equilibrium pa- 
rameters of isolated planar domain walls and derive the 
equilibrium parameters of stripe domains in system with 
out-of-plane magnetization. 



A. The structure of domain walls 

For a planar isolated domain wall with energy density 
3>(0) Q the structure is derived by optimization of the 
functional 



Xo 



WDw = Adl + K^A^{9) 



(17) 



with the boundary conditions 9x{ioo) = 0, 9{+oo) = 6i, 
9(—oo) = 92 {x is a spatial variable across the domain 
wall), and A$(6') = $(6*) - $o, where $o = ^(^i) = 
'^{92) = {min<f>} is the global minimum of the system. 
For such a one-dimensional problem ( 17 1 the Euler equa- 



tion and the first integral can be written as [43] 

2xl9xx = 
xl{9xf = A^{9) 



(18a) 
(18b) 



where xq ~ \/ AjK^ is a charactristic width of the do- 
main wall. 

The domain wall profiles 9(x), their energy and charac- 
teristic sizes can be readily derived by direct integration. 
However, Eqs. (18a I allow to understand the main fea- 
tures of such solutions. The Eq. (18a I shows that the 



inflection points of the domain wall profiles 9{x) corre- 
spond to stationa ry po ints of potential ^ {^g — 0). The 
second equation (18b) shows that the larger the devia- 



tions of the energy from the minima A<i>(6') the larger 
the gradients of the profiles, 9x- Thus, the nucleation 
and further evolution of local minima in the metastable 
region causes complex reconstructions of the domain wall 
profiles. Transformations of the domain wall profiles have 
been earlier observed in easy-axis antiferromagnets and 
other magnetic crystals (see examples in Ref. Due 
to the remarkable lability of the potential profiles (|4| 
(Figj2]) this effect is expected to be strong in the systems 
with competing anisotropies. As an example, we derive 
the parameters of the domain walls for four-phase do- 
mains with the canted states ( 15 ) that are realized in the 
systems with a = tt/A, l>>!:>Oat zero fields. 

In this case three types of domain walls can exists 
(Fig.[8];A)): DWI between domains with 6*1, 6*2 and 613, 6*4 
{A9i = - 6*2! = |6'3 - 6*4! = arccos>f), DW II between 
9i, 93 and 92, 9^ {A9u = l^i - ^sl = 1^2 - = tt- A^i), 
and DWIII of 180 type between domains 9i, 6*4 and 6*2, 
By integration of ( |18b ) the energy a and the mag- 



netization profiles for the DWI (upper signs) and DWII 
(lower signs) can be readily obtained as 



^2(1 - >^') 



In T 



tan 9 — tan ( 
tan 9 + tan ( 



(20) 



So \/l — T ^ arccos ( ± xr) 



where 6q = ^/AKc- Depending on the ratio 1 > x: > 
DWI with A9i < 90 becomes more favourable than 
DWII with A9jj > 90 and should exhibit stronger con- 
trast in experiment (71!. Such domain boundaries were 
experimentally observed in thin films of (Ga,Mn)As by 
Lorentz microscopy [71] . From Refl72l and using ( pO| ) we 
can evaluate the DW width. For DWI at the temperature 
T = 10 (30)K we obtain S = 50 (lOO)nm while for DWII 
(5 = 43 (76)nm which is consistent with the experiment. 
Here, A = 0.4 • lO^^^Jm"!, if^ = 1-18 (0.32)Jm-3^ Ku = 
0.18 (0.11)Jm-3. For x > 1 only DWIII exist. 

These domain walls are characterized by strong varia- 
tion of their parameters with the applied magnetic field. 
At magnetic field or domain walls of two types 
are present: walls where the magnetization vector ro- 
tates less or more than 180° (Figj8];A),(B)). Note, that 
for magnetic field > 2(1 — x) the metastable min- 
ima (Figj2jB)) strongly influence the profile and energy 
of the domain. Such a remarkable modification of the 
structure should strongly infiuence magnetoresistence of 
domain walls (e. g. in nanoconstrictions) . 

In Fig|8jC),(D) domain wall profiles and typical phase 
portraits for the case a = (xr = 0.7) are plotted. At 
applied magnetic field only 360 domain walls exist. 
Within these walls, nuclei of the domain with tt and ±7r/2 
are present (see phase portrait in Figjsj^C)). In a magnetic 
field strong enough these nuclei disappear. So, the energy 
of domain wall increases, although the width decreases. 
By application of magnetic fields perpendicular to the 
easy axis a in the interval [0; oi] there are walls of two 
types between upper and lower canted phases (Fig|8|D)). 
In each domain wall nuclei of domains with 7r/2 and 37r/2, 
correspondingly, are present. At ~ hj.{a,i) these nuclei 
expand forming three-phase multidomain textures |27j . 



B. Parameters of stripe domains 

Magnetic configurations in (Ga,Mn)As systems include 
a number of noncoUinear two- and multi-phase states. 
These phases can create thermodynamically stable mul- 
tidomain states |42l . For two coexisting phases with the 
magnetization M^-'^' and M^^^ effective values of magne- 
tization can be introduced |42] 

Af_L = (M(i' -M(2)) .v/2, 

H= [h - Ht^ - 47r(M(i) + M*^))] • v/2. (21) 

In particular, for perpendicular magnetized (Ga,Mn)As 
nanolayers with k > b the phase diagrams of magnetic 
states (Fig. [ll(A),(E)) are similar to those for uniax- 
ial ferromagnets. In this case M±_ = Mq and domains 
are separated by 180° domain walls [43]. According to 
Eq.( 21 1 the problem of multidomain states for two-phase 



(19) 

noncoUinear states can be reduced to a ferromagnetic 
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FIG. 8; Domain wall profiles for x = 0.7: (A), (B)- a — n/A, (C), (D) - q = 0. Corresponding energy profiles are plotted in 
Fig. [2] Insets show the phase plane (6; Ox) where 9^ ~ dO/dx (See Eq.(18l). 
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FIG. 9: Characteristic lengths li(ii) for stripe domains with 
DWI(DWII) in four-phase state (Zo = -R-ZAKl/iAMl). 



coUinear domain structure with up and down magneti- 
zation Mj_ in a bias field H (v is the unity vector per- 
pendicular to the layer plane, Ht^ is the transition field 
between the phases M^^^and M^^^) Similar mul- 

tidomain texture are formed in the systems with a num- 
ber of coexisting phases larger than two. For example, 
for —\<x<\,a — Q the magnetic configurations 
(Eq ~ 



(15)) create four-phase multidomain states. In 
a layer with a||v these textures can be described by a 
model of ferromagnetic domains with the magnetization 
Afj^ = M^{1 + k)/2. (Fig. Igl. 



With effective values of the magnetization Mj_ and bias 
field H (21) the energy density of a (Ga,Mn)As nanolayer 
with stripe domains can be reduced to the well-studied 
model of ferromagnetic stripes [331 [731 [23] 



w 



2'KMi 



w„ 



2Ap 



Hq 



2nM 



± 



with the stray field energy density given by 



w„ 



2p 



•7r2 , 
^ -10 



(1 - t) In 



cos^ {'Kq/2) 
sinh^ {pt/2) 



(22) 



dr , (23) 



where p — 2Tit/D, q = {d+ — d^)/D, t is the layer of 
thickness t, D = d^ + d- is the stripe period, and d± are 
domain sizes with up and down magnetization (Fig. |9|. 
The dimensionless parameter A = a{>c) / {4:M±{K)'^t) = 
T:l{x))/t measures the ratio between the domain wall en- 
ergy (7 and the stray field energies. It scales with the 
characteristic length 1{k) — a/{4iTM'^) as the relevant 
material parameter 23] • Minimization of (22) with re- 
spect to p and q derives the solutions for the geometric 
parameters d± as functions of three control parameters 
of model (23): the layer thickness t, the bias field H, and 



factor A. These solution have been investigated in detail 
(see (431 [731 [74] and bibhography in Ref. [^). Partic- 
ularly, it was shown that the solutions for stripes exist 
only below certain critical field H < H*{A) < ATrM± 
[74] . As the bias field approaches H* the stripes gradu- 
ally transform into the homogeneous state by unlimited 
growth of the period (D 00). However, at the critical 
field the domain of the minority phase preserves a finite 
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size d-{H*) = d*_. At higher fields {H* > H > 47rMj_) 
it exists as a metastable state gradually shrinking to zero 
size at the saturation field. 

In perpendicular magnetized (Ga,Mn)As layers the pe- 
riod of multidomain patterns exceed their thicknesses 
[5U1 [551 155] , For such large stripes {D > t) the expan- 
sion of the integral ( p3| allows to simplify the problem 
[73] , After some algebra, the solutions for stripes can be 
derived in analytical form as a set of parametrical equa- 
tions 



D{H) = ■Kut/\Jl - {H/H*f, 
d± = (D/t:) HTCCOsiT H/H*), 



(24) 



H*{u) = 4TTMj_f{u), 2A = g{u). (25) 
Here we introduce parameter u — d*_/t and functions 

f{u) = [2arctanl/'u - uln(l + 1/u^)] /tt, 

g{u) = {l + u^) ln(l + u^) ~ V? (26) 



According to (24 1 at zero field -D(O) ~ ~ ntu, thus, 
the ratio Do/d*^ ~ tt. It means that at the transition 
field the domain of the minority phase becomes approx- 
imately six times narrower than the domain size at zero 
field {Dq/2). Within this approximation the equilibrium 
magnetization in the stripe phase equals 



(Af) Mj_q = (2A/_L/7r) arccos(iJ/i?*). 



(27) 



Finally, in the limit of large domains D ^ t, Dq = 
7rfexp(A — 1/2), and the transition fields for stripe and 
bubble domains becomes exponentially small, e. g. tran- 
sition field H* = 4M_L exp(-A-|-l/2), the bubble collapse 
field Ha = 16Mj_exp(-A - 1/2), and ratio H*/Ha = 
e/4 = 0.6796 [73]. 

These results show that the solutions of magnetic 
domains in (Ga,Mn)As layers should demonstrate gen- 
eral features similar to those in uniaxial ferromagnets. 
However, there is an important difference between these 
two systems. In uniaxial ferromagnets the characteris- 
tic length is expressed as a combination of basic mag- 
netic parameters (constants of uniaxial anisotropy K, 
exchange stiffness A and saturation magnetization M): 
If = cr//(47rM2) = y/AK/{nA'P) depends only on val- 



^ . . . 

ues of uniaxial anisotropy K. On the contrary, in the di- 
luted magnetic semiconductors the characterstic length 
l(>c) strongly depends on the values of competing mag- 
netic anisotropics and varies in a broad range providing 
a complex behaviour of multidomain patterns in these 
materials. Eqs. p4|, (p5|, (27 1 connect equilibrium pa- 



rameters of stripe domains with material parameters of 
(Ga,Mn)As systems Q. For this model calculations of 
M± and the domain wall energy a allow to express A 
as a function of k and a. For example, for four-phase 
domains (Eq. (|15[)) the magnetization and the domain 



wall energies (Eq. (19l) yield functions Z/(//)(>f) plotted 
in Fig. |9] 

Stripe domains have been observed in a number of 
(Ga,Mn)As nanolayers [SSEJEIIMIIMIIMIISSII^- In 

Table I we collect experimental data (indicated by bold) 
and the calculated stripe domains parameters (by solv- 
ing Eqs. (p4l), (p5l), (pTl)) for (Ga,Mn)As layers (a-c), an 



yttrium-iron garnet film (d) and FePd nanolayers (e,f). 
For thin layers of (Gao.957Mno.o43)As for T = 9 K (1) 
and T = 80 K (2) [SU] we use the experimental values of 
t and -Do to calculate other parameters of stripes. For 
a layer (Gao.93Mno.o7)As at T = 80 K [55] we use t and 
a value of the transition field H* /{4:TtM±) to calculate 
I, A and domain sizes Dq and c?^. For comparison we 
derive stripe domain parameters for an epitaxial garnet 
film Yi.88Luo.2Cao.92Geo.92Fe4.o80i2 [7Sj, and for FePd 
nanolayers [75] . According to [55] the saturation mag- 
netization in (Gao.93Mno.07) As Ms = 28 kA/m. Then, 
from the calculated value of the characteristic length we 
derive a — 49.0 /iJ/m^ (for comparison, in the (Y,Fe) 
garnet film = 13.6 kA/m, and a = 110 ^J/m^ [75j). 

In garnet films and other classical materials with per- 
pendicular anisotropy regular stripe domains are ob- 
served, if the layer thickness is considerably larger than 
the characterisitic length. In such systems the equilib- 
rium domain sizes at zero field Do do not exceed the 
layer thicknesses {Dq < t). When films becames thin- 
ner than / (e.g. in the vicinity of the compensation 
temperature of ferrimagnets [73 ) the demagnetization 
forces are too weak to overcome coercivity and the for- 
mation equilibrium domains is impeded. As a result 
such multidomain patterns consist of very large domains 
with irregular boundaries |43j . Similar disordered do- 
mains and strongly hysteretic behaviour have been ob- 
served in (Ga,Mn)As films with perpendicular anisotropy 
[SnilSl IMJ Ell IM]- It means that the equilibium mul- 
tidomain state are hardly reached in these systems. For 
example, in (Gao.93Mno.o7)As layers the observed width 
of the minority stripe at the critical field is d*_ = 1.7 ±0.3 
/im '66], and about one order larger than the calcu- 
lated equilibrium value (dl = 0.2 /zm). Up to now only 
few results on experimental investigations of multido- 
main states in perpendicularly magnetized (Ga,Mn)As 
nanolayers have been reported. New detailed investiga- 
tions involving modern experimental methods developed 
in other fields of nanomagnetism (see e.g. [TT]) are desir- 
able. The results of this section 

establish important physical connections with multido- 
main states in other classes of perpendicular magnetized 
materials and provide a theoretical basis for future re- 
search. 



VII. SUMMARY AND CONCLUSIONS 

We have developed some micromagnetic methods [42l 
I43| which give a consistent description of magnetization 
processes and multidomain structure in systems with 
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TABLE I: Parameters of the stripe domains in (Ga,Mn)As layers (a-c), a (Y,Fe) garnet film (d) and FePd nanolayers (e,f). They 
include experimental values (bold) and the results derived from model (22 I. Here, t is the layer thickness, / is the characteristic 



length, A = nl/t is the dimensionless parameter measuring the ratio between the domain wall energy and the stray field energy 
(see Eqs. (22), (23))), Do is the equilibrium period at zero field, dl is the equilibrium size of the minority phase a.t H = H* , 
the trasition field into the homogeneous state, Eq. (25). 







t, fim 


I, fiin 


A 


Do, ^m 




J/7(47rMx) 


(a). 


(Gao.g57Mno.o43)As |30 


0.2 


0.132 


2.07 


3.0 


0.95 


6.6 10-" 


(b). 


(Gao.957Mno.o43)As [30 


0.2 


0.220 


3.45 


12.0 


3.82 


1.7 10"" 


(c). 


(Gao.93Mn().o7)As 66^ 


5 10"" 


0.10 


1.920 


0.643 


0.2038 


3.9 10 " 


(d). 


(Y,Fe) garnet film ES] 


11.0 


0.47 


0.1342 


12.50 


3.0 


0.6 


(e) 


. FePd nanolayer [76] 


3.6 10 " 


9 10"^ 


0.7526 


0.13 


3.9 10~" 


0.2612 


(f) 


. FePd nanolayer [7B] 


1.15 10"" 


9 10"^ 


2.477 


0.26 


8.3 10"" 


4.4 10"" 



competing anisotropies such as diluted magnetic semi- 
conductors. 

Theoretically constructed phase diagrams in external 
field components in the limiting case of ideally soft mag- 
netics allow to understand the creation of equilibrium do- 
main structure |29j . within ahysteretic magnetization re- 
versal, and explain various parts of magnetization curves 
(Figj?]). Thus, magnetic phase diagrams allow to put 
in good order and classify a vast amount of experimental 
data on reorientation effects, multidomain processes, and 
magnetization reversal in (Ga, Mn) As systems. These 
diagrams also give opportunity to predict change of mag- 
netic states of the system in zero magnetic field under 
action of temperature [301 1311 132 and in the case of arbi- 
trary angle a between competing anisotropy axes |24l I25j . 

It is also shown that the applied magnetic field causes 
drastic transformations of the domain wall profile and 
strongly infiuences its parameters. Domain walls can 
serve as nuclei of domains for a new phase. At certain 



values of the magnetic field a domain wall can be divided 
into domains of a new phase and two types of new domain 
walls. At certain critical endpoints of phase coexistence, 
domain walls can disappear by the rotation of the mag- 
netization in adjacent domains towards each other. 

For nanolayers with perpendicular anisotropy the geo- 
metrical parameters of stripe domains have been calcu- 
lated as functions of a bias field. 
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